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Abstract — In  the  dynamic  electromagnetic  environments,  a 
polarization  estimator  is  usually  utilized  by  an  adaptive 
polarization  filter  to  acquire  the  knowledge  of  the  polarization  of 
the  interference  field,  then  optimal  polarization  is  calculated  for 
polarization  filtering.  In  this  sense,  an  adaptive  polarization 
filter  can  be  modeled  as  an  open  loop  cascading  system  of  a 
polarization  estimator,  an  optimum  polarization  calculator  and 
a  polarization  filter.  Based  on  the  dual  channel  orthogonal 
polarization  measurements,  theoretical  performance  analysis  of 
the  polarization  estimator  is  derived.  It  is  proved  that  the 
polarization  estimation  error  obeys  the  Rayleigh  distribution 
approximately.  Finally,  by  studying  the  transmission  process  of 
the  estimation  error  in  a  polarization  filter,  a  series  of  analytic 
formulae  on  the  performances  of  adaptive  polarization  filters 
are  derived,  which  statistically  reveal  the  relationship  between 
the  filtering  performance  and  the  signal  to  noise  ratio  clearly. 

Index  Terms — radar,  polarization,  estimation,  polarization 
filter,  adaptive 

I.  INTRODUCTION 

Under  the  condition  of  modern  electronic  warfare, 
environments  become  more  and  more  adverse  for  the  radar. 
The  problem  of  anti-jamming  is  of  great  importance.  The 
adaptive  polarization  filtering  technique  has  been  considered 
as  a  promising  technology  for  radar  anti-jamming  [1-3,8, 9]. 
In  the  past  thirty  years,  lots  of  adaptive  polarization  filters 
have  been  designed  for  various  purposes.  Many  of  them 
worked  well.  However,  their  intrinsic  limitations  were  also 
exposed.  Hence  evaluating  the  performance  of  an  adaptive 
polarization  filter  becomes  urgent  for  the  polarization  signal 
processing  technique.  Due  to  the  complexity  of  clutter,  the 
diversity  of  interference  and  the  specialty  of  polarization 
filtering,  the  designation  of  polarization  filters  usually  has 
very  strong  pertinence.  There  is  still  no  frame  theory  on  the 
performance  of  polarization  filtering. 

Essentially,  polarization  filters  improve  the  qualities  of  the 
desired  signals  through  selecting  the  antenna’s  polarization. 
There  are  mainly  two  classes  of  polarization  filters.  One  is 
called  Interference  Suppressed  Polarization  Filters(ISPF) 


which  makes  the  received  jamming  power  minimum 
through  selecting  antenna’s  polarization  orthogonal  to  that  of 
jamming.  No  minding  the  signals’  polarization  is  one 
characteristic  of  these  filters.  Minimizing  the  jamming  power 
is  the  only  goal  of  ISPF [1,2,9].  The  other  class  filters  such  as 
SINR  (Signal  to  Interference  plus  Noise  Ratio)  and  PDSI 
(Power  Difference  of  Signal  and  Interference)  filters  consider 
the  polarization  of  both  signals  and  jamming  in  the  antenna’s 
beam[4’7]. 

In  order  to  improve  the  quality  of  the  received  signals,  the 
prior  information  about  jamming  must  be  known.  In  fact,  the 
polarization  of  jamming  is  usually  obtained  by  estimation 
[4,6].  As  for  jamming  of  high  polarization  degree,  jamming 
suppression  ratio  can  be  very  large  (infinitude  in  theory). 
However  if  the  receiving  polarization  isn’t  orthogonal  to  that 
of  interference  strictly,  the  filtering  performance  will 
decrease  rapidly  [2].  Hence,  the  estimation  to  interference 
polarization  is  a  key  in  the  filtering  process. 

The  physical  premise  of  obtaining  the  polarization  of 
incident  electromagnetic  field  is  the  dual  channel  orthogonal 
polarization  measurements.  Theoretically,  the  amplitude  and 
phase  of  the  measurements  can  represent  the  polarization  of 
incidence  field  uniquely  [8].  In  practice,  the  outputs  of  the 
dual  channel  orthogonal  polarization  measurements  are 
taken  as  the  incident  waves’  estimations,  that  is  to  say,  the 
dual  channel  orthogonal  polarization  measurements  can  be 
regarded  as  a  polarization  estimator.  From  the  viewpoint  of 
signal  estimation,  it  is  an  optimum  polarization  estimator, 
and  it  represents  a  superior  limit  of  filtering  performance  [13]. 

This  paper  is  organized  as  follows.  In  Section  II,  the 
model  of  the  polarization  estimator  is  established.  Its  error 
based  on  the  Stokes  vector  is  derived.  In  Section  III, it  is 
proved  that  the  polarization  estimation  error  approximately 
obeys  the  Rayleigh  distribution.  In  Section  IV,  a  series  of 
analytic  formulae  on  the  performances  of  adaptive 
polarization  filters  are  derived,  which  stochastically  reveal 
the  relationship  between  the  filtering  performance  and  the 
signal  to  noise  ratio  clearly. 
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II.  The  Polarization  Vector  Estimation  Based  on 
the  Dual  channel  Orthogonal  Polarization 
Measurements 


Select  horizontal  and  vertical  linear  polarization  (h,  v) 
as  polarization  basis,  and  assume  any  a  planar  incident  wave 
can  be  represented  as  S  =  [ sh ,  sv  ]  based  on  ( h ,  v) .  Here 

the  symbol  “  T”  represents  transpose.  Suppose  this  wave  is 
received  by  dual  channel  orthogonal  polarization 
measurements.  Then  the  output  is  a  complex  vector  and  is 
written  as  X  .  Neglecting  the  dual  channel  inconsistency, 
coupling  and  nonlinearity,  the  signal  model  can  be 
represented: 

x  =  s  +  n  (1) 

Where  n  is  the  thermal  noise  of  the  measurement  system.  In 

fact,  the  output  is  directly  looked  as  the  estimation  of  the 

incident  signal  vector,  that  is: 

*  x  s  +  n 

s  = - = - 

x  ||  ||  s  +  n 

Where  S  denotes  the  estimation  of  the  incident  wave,  1 1  •  1 1 
is  the  2  norm  of  complex  vector.  Obviously  ||  S  ||=  1 .  To 
simple,  suppose  the  power  of  the  signal  is  unit,  that  is 
1 1  S  1 1=  1 .  Hence  S  can  be  written 

s=[cos7,  sin  ye*J,ye[Q|],0e[Q2rc]  (3) 

Where  (7, 0)  is  called  polarization  phase  descriptor [8]. 

In  order  to  analyze  the  performance  of  the  estimator,  the 
difference  between  the  estimation  polarization  s  and  the  real 
polarization  s  must  be  measured.  From  the  relationship  of 
the  Stokes  vector  and  electric  field  vector,  it  is  known: 

Js  =  R(s®s *) 

Js  =  R(s®s *) 

Where  J  s  and  J  s  are  Stokes  vector  of  the  true  polarization 

and  estimation  polarization,  respectively.  The  symbol  “  ®  ” 
and  “  *  ”  denote  Kronecker  product  and  conjugating 
respectively  .  R  is  written  as:  [8]. 

"lOO  1  " 

10  0-1 

R  = 

0  110 
_0  j  -j  0  _ 

Since  ||  5  I  HI  S  ||=  1,  Js  =  ru;  f  Js  =  [l,gTs  ]T  ,  in 
which  gs  and  gs  are  called  the  Stokes  Subvector  of  true 
polarization  and  estimated  polarization.  Obviously,  gs  and 
gs  are  both  real  three  dimension  vector.  In  addition,  they 


satisfy  ||  gs  |  |=|  |  gs  |  |=  1  and  are  mapped  to  two  points  on 
the  Poincare  polarization  sphere[8],  as  illustrated  in  Fig.l. 


Fig.  1  the  vector  presentation  of  true  polarization  and  estimated  polarization  on  Poincare  sphere 

Suppose  6  is  the  angle  between  true  polarization 
gs  and  estimated  polarization  gs  .  It  is  the  natural  scale  for 

the  performance  of  polarization  estimator.  It  can  been  easily 
proved  that: 

2(1  +  cost?)  —\  sH s  \2  (4) 


The  symbol  “  H”  is  complex  conjugate  transposition.  It 
shows  that  the  angle  between  true  polarization  and  estimated 
polarization  can  be  calculated  by  the  inner  product  of  electric 
field  vector.  Substitute  Equ.  (2)  into  above  formula,  then  we 
get: 


sm 


v2y 


n 


nHs  I2 


(5) 


1  +  n  s  +  s  ii+  ||  n  || 

Suppose  the  SNR  of  measure  system  is  high,  that  is 
||  il  ||«||  S  ||—  1 .  Then  Equ.  (5)  can  be  rewritten  as: 


sm 


nHWn 


(6) 


Where  W  =  I2  —  ssH  ,  /2  =  diag{ I,l}.  Since  ||  s  ||=  1 , 
the  matrix  W  must  be  nonnegative  definite  Hermitian  matrix. 
In  fact,  substitute  Equ.  (3)  into  the  expression  of  W,  then 
Unitary  diagonalize  W,  we  get: 


w  =  uhau 


(7) 


in  which  A  =  diag{  1,0}  ,  U  is  a  Unitary  matrix,  its 


expression  is 

u=  sin  jej*  -cos  7  (8) 

cos  7  sin  ye~ji * 

r\ 

Substitute  Equ.(2)  into  Equ.(6),  then  sin  2  —  ~  m  H  Km  • 

2 

Where  m  =  Ufl  is  the  Unitary  transform  of  the  noise  output 
vector.  Since  ||  ft  ||=||  ftl  ||«  1 ,  then  6  «  1 ,  then 

0  ~  2^1  mH Am  (9) 

From  Equ.(6)  and  Equ.(9),  the  6  is  dependent  on  two 
factors:  one  is  the  system  output  noise,  another  one  is  the 
signal  polarization  parameter  whose  effect  is  embodied  in 
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matrix  W  According  to  Equ.(9),  we  will  analyze  the 
statistical  property  of  polarization  estimation  error. 


III.  Performance  Analysis  of  the  Polarization 
Estimator 

For  radars  operating  above  VHF/UHF  frequencies, 
regardless  of  clutter  and  interference,  the  output  noise  vector 
n  of  dual  channel  is  mainly  resulted  from  reception  thermal 
noise  which  can  often  be  viewed  as  band  limited  white  noise 
and  obeys  Gauss  distribution.  Suppose  tl  is  a  complex 
Gauss  vector  with  zero  mean  and  covariance  matrix 

—lnnH^j  ,  Here  the  symbol  of  denotes  statistical 

expectation.  Denote  m  =  Ufl  ,  so  m  is  a  complex  Gauss 
vector  with  zero  mean  too  and  its  covariance  matrix  is 
£  m  =  UZnU  H  - 

Denote  m  =  [mx,m2]T  =  [xx  +  jyl,x2  +  jy2Y , 
mi  is  a  gaussian  complex  variable,  according  to  stochastic 
process  theory,  mi  is  a  narrowband  gaussian  noise  actually,  x\ 
andyi  are  two  real  low  frequency  orthogonal  process [11].  If 
the  covariance  of  m\  is  GmU  ,  then  the  covariance  ofxi  andyi 


are  <Jmll/2  .  The  envolope  of  m\  Ax  =  xf  + 
satisfies  Rayleigh  distribution,  and  its  PDF  is 


yi 


2x 

A  (*)  = - exP 


.  2  ') 


X  >  0 


Substitute  m  into  Equ.(9),  then  0  ~  2  Ax .  So  0  obeys  the 
Rayleigh  distribution  approximately  and  its  PDF 


f(0)  = 


If 


9 


2(7 , 


-exp(  - 


9 


4(7, 


-),  9  e  [0 ,n] 


V,  = 


(10) 


(11) 


According  to  the  relationship  between  Zm  and  , 

=d,n  sin2  r+^,22  cos2  7_2Re[c7(il2e/^]sinycosy  ( 1 2) 

Especially,  when  dual  channel  are  independent  to  each  other, 
that  is  to  say,  the  two  polarization  components  of  n  are 
independent  to  each  other,  then  an]2  =<7*21  =0 .  Based  on 
Equ.(12),  thencr7Mll  can  be  written  as: 

<7*11  =^„iisin2r  +  0'B22cos27 
Furthermore,  if  the  power  of  the  two  polarization 
components  of  it  are  equal,  namely,  Gnl  x  —  On22  ,  then 


GmX  |  can  be  rewritten  as:  GmX ,  =  GnX  x .  Fig.2  offered  the  PDF 

curves  of  6  based  on  two  different  SNR.  In  the  figure,  SNR 
is  the  Signal  to  Noise  Ratio  in  the  output  of  dual  channel.  The 
stick  curve  is  statistical  graph  of  PDF  of  from  20000  Monte 


Carlo  experiments.  Solid  line  is  the  theory  curve  based  on 
Equ.(10).  From  Monte  Carlo  experiments,  in  higher  SNR  ( 
>=10dB  )  ,  the  Rayleigh  distribution  can  be  good  in 
representation  of  the  PDF  of  6  . 


Fig  2  the  PDF  curve  of  polarization  estimation  error  angle  6 
(The  signal  polarization  is  Left  Circularly  and  the  two  outputs  of  dual  channel  are 
independent  and  their  powers  are  equal) 


IV.The  Performance  Analysis  of  Adaptive 
Polarization  Filter  Based  on  Polarization 
Estimation 

The  technique  of  polarization  filtering  must  know  the 
polarization  of  interference  filed  beforehand.  Actually,  many 
polarization  filters  in  applications  have  an  adaptive  estimator, 
such  as  Nathanson’s  APC  cancelled 14],  Gherardelli’s  MFP 
APC  and  MFP  SAPC  canceller[15,16] ,  Pother’s  optimum 
polarization  detector[5]  and  so  on.  These  filtering  often  have 
three  steps:  at  first,  adaptive  estimate  the  interference 
polarization,  then  calculate  optimum  receiving  polarization, 
finally,  filter  with  interference  suppressed.  Obviously,  the 
estimation  precision  of  interference  will  directly  affect  the 
performance  of  the  whole  filtering  process. 

The  characteristic  of  polarization  estimator  error  is 
analyzed  in  above  two  Sections,  its  influence  on  the 
performance  will  be  studied  in  this  Section. 

Denote  the  estimation  polarization  is  S  .  gs  and  gs 
are  respectively  Stokes  Subvectors  of  true  and  estimation 
polarization,  and  the  angle  between  the  two  vectors  is  6 

which  satisfies  =  .  If  there  is  only  one 


interference  source,  then  the  process  of  polarization  filtering 


can  be  represented  as  in  Fig  3. 


Incident 

wave 


/ave  ^ 


Estimation  interference 
field  polarization  s 
based  on  dual-channel 


Figure  3  the  model  of  open  loop  cascading  system  of  a  adaptive  polarization  filter 
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Define  s  =  [cos/,  sinfeJ^  J  ,  then  its  orthogonal 

polarization  is  S±  =  [—  sin  f,  COS  fe  ^  ]  .  Obviously,  its 
Stokes  Sub  vector  satisfies  Si  =  ^nSs  .  Where 
A12  =  diag\—  \—  l,l}  .  For  the  limit  of  polarization 
estimation  precision,  s±  isn’t  always  orthogonal  with  S 

rigidly.  That  is  to  say,  there  must  be  a  little  interference  signal 
leak  into  the  radar  receiver,  it  is  called  interference 
remainder: 

Pr  =1*1*  l2=  \jlu  *Js  =  y(l  + 

Where  ,7  ~r{s±  [l,  gTL  J  is  Stokes  vector  of 

S±  correspondingly,  U4  =dia$X\~ty  ,  A3  =  diag{  1,1,— 1}  . 


Because  of  g±  =  Augs  and  COS0  =  g^ gs  then  formula 
above  can  be  written  as: 

p,  =  |  [l +(Ala#,  r  ]=  |  (l  -grs, ) = 1  (1  -cose)  <  13) 

It  is  a  very  significant  conclusion.  It  indicates  that 
interference  remainder  is  only  decided  by  the  estimation  error 
of  polarization  for  all  adaptive  ISPF.  It  also  proves  that  the 
polarization  estimation  is  the  key  of  polarization  filter  design. 
In  addition,  it  presents  a  simple  mathematical  model  for 
adaptive  polarization  filter  performance  analysis.  For 
example,  we  can  easily  get  the  probability  property  of  Pr 

only  if  the  statistical  property  of  9  is  known.  In  the 
following,  some  analytic  formulae  will  be  deduced.  They  are 
universal  for  all  adaptive  polarization  filters  and  offered  a 
super  limit  for  their  performance. 

When  the  interference  is  strong  and  thermal  noise  is  weak, 
polarization  estimation  error  angle  6  obeys  Rayleigh 
distribution  approximately,  and  its  PDF  is  offered  by  formula 
(10).  Because  the  relationship  between  interference 
remainder  Pr  and  0  is  monotonously  descending  function, 


then  the  PDF  of  P„  can  be  derived  from  formula  below. 


fr,(n)  =  fo(0) 


de 


drj 


=  fe(0)~ 


sin  6 


where  Tj  =  Pr  .  Substitute  6  =  arccos(  1  —  2 77 )  and 

sin  0  =  2  ^77  -  77 2  into  above  formula,  then 

r  ,  arccos(l  -  2?7 )  [  [arcco s(l  -  2r\  )]2  ) 

fniri)m - v  /  exp  ^  - '-U-\,TJ  e  [0,  1] 

2  amliy/ri -ri2  {  4crmll  J 

(14) 


suppression  ratio”.  For  example,  a  narrowband  Left  Handed 
Circularly  polarization  interference  signal,  the  dual  channel 
estimator  is  adopted  in  radar  processing  and  polarization 

estimation  error  angle  is  6  =  5° .  Based  of  the  estimator,  use 
orthogonal  estimation  polarization  to  suppress  the 
interference,  according  to  Equ.(13),  we  obtains  7]  —  0.002  . 
Namely,  the  polarization  interference  suppression  ratio  is 
-27.2dB  .  In  some  sense,  Equ.  (14)  is  the  PDF  function  of 
polarization  interference  suppression  ratio  Tj  of  polarization 
filter.  Fig.  4  presents  some  PDF  curve  of  Polarization 
Interfere  Suppression  Ratio.  It  indicates  that  the  higher  the 
SNR,  the  more  7]  approaches  zero,  in  other  words,  the  better 

the  performance  of  the  filter  is. 


Fig.  4  the  PDF  curve  of  Interference  Suppression  Ratio 

From  the  view  of  statistics,  the  PDF  of  polarization 
interference  suppression  ratio  is  a  whole  representation  for 
the  performance  of  interference  suppression  polarization 
filter.  But  in  practice,  using  one  or  more  statistical  scale 
representing  the  performance  is  more  popular.  Similarity  to 
the  concepts  of  “Trust  Level”  and  “Trust  Interval”  in 
statistics[17],  we  choose  “probability  level”  representing  the 
performance  of  interference  suppression  polarization  filter: 
given  a  probability  level  OC ,  it  corresponds  to  a  interference 

suppression  ratio  critical  value  7]0  ,  interference  suppression 
ratio  of  polarization  filter  7]  should  be  smaller  than  77 0  in 
probability  of  OC ,  that  is 

JJ°  /„  (l)dn  <  a  (15) 

Substitute  Equ. (14)  into  above,  then  we  get: 

no  =  j  t  - cos  V-  4°™i  1  ln(t _  a)  1  06 

) 


Since  1 1  S  \  |=  1 ,  the  remainder  power  of  interference  Pr 
( or  T]  )  must  be  between  0  and  1 ,  hence  Tj  can  be  regarded 

as  the  attenuation  of  interference  signal  through  the 
polarization  filter,  or  be  called  “polarization  interference 


Fig  5  presented  the  curve  of  critical  polarization 
interference  suppression  ratio  v.s.  SNR  and  v.s.  probability 
level.  Fig  5(a)  indicates  that  the  relationship  between  Tj0  and 
SNR  is  mono-descending  with  condition  of  a  given  a 
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probability  level  OC .  And  when  SNR  is  higher  than  5dB,  the 
relationship  between  7]0  and  SNR(dB)  is  linear 
approximately.  Fig  5(b)  indicates  that  the  relationship 
between  7]0  and  probability  level  OC  is  mono-increasing  on 
condition  of  a  given  SNR. 


Figure  5  the  curve  of  Critical  Polarization  Interference  Suppression  Ratio  v.s.  SNR  (a) 
or  probability  level  (b)  :  (a)  From  up  to  down,  probability  levels  respectively  are  OC  =0.99 
,  0.9,  0.8,  0.7,  0.6,  0.5  (b)  From  up  to  down,  SNRs  respectively  are  SNR=  1 0 ,  15,  20 
,  25dB 


Define  SNR  =  5 


If  )  is  the  signal  to  noise 


ratio  of  actual  output  of  dual  channel,  where  the  bottom  label 
“g”  denotes  the  actual  output  of  the  measurement  system. 
From  the  suppose  that  the  interference  field  has  unit  power, 

"  ^  —  l/ truL  /||s||  ),  where  the 

n’Z  /  'n’g  II  II je  / 


then  SNR  =  5 


tfLn 

g  /  n’g 


symbol  of  “  tr”  represents  the  trace  of  matrix[10].  It  can  be 
concluded  that  the  covariance  matrix  £  of  n  is  actually  the 


standardization  of  the  covariance  matrix  g  to 
interference  power,  that  is  =  Zwgy/||s||^  .  Similarly, 

’  ^*11  ■  Especially, 

when  the  noise  of  dual-channel  have  the  same  characteristic 
and  in  independent,  then£w  g  =  G  nXX  gI2  •  Since  it  is  known 

1 

before  GmXXg  =  OnXXg  =  -trl,ng,  then 
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1 

2  SNR 


(17) 


Substitute  Equ.(17)  into  Equ.(10),  (14)  and  (16),  then  some 
relationships,  which  include  the  PDF  of  Polarization 
estimation  error  angle  9  and  polarization  suppression 
ratio  77  ,  and  the  function  between  Polarization  Interfere 


Suppression  Ratio  critical  value  7 ]0  and  Signal  to  Noise  Ratio 
(SNR),  are  all  derived. 


f0(6)  =  SNR  •  9  exp 


f  SNR  02^ 


,  de[0,n] 

(18) 


L(n) 


SNR  •  arccos(l  -  2r])  f  SNR 


yfn- 


R 


exp<j - —  [arccos(l  -2 rj)]2\,  r/e  [0,1] 


(19) 


From  these  formulae,  the  relationship  between  the 
performance  of  estimation  and  Signal  to  Noise  Ratio  can  be 
clearly  known. 

V.  Conclusion 

An  adaptive  polarization  filter  can  be  modeled  as  an 
open  loop  cascading  system  of  a  polarization  estimator,  an 
optimum  polarization  calculator  and  a  polarization  filter. 
Based  on  this  model,  a  super  limit  of  theoretical  performance 
of  the  adaptive  polarization  filters  was  obtained  and  its 
analytical  formula  is  presented  too.  The  study  is  of  great 
importance  to  the  performance  elevation  of  adaptive 
polarization  system  analysis  and  optimum  design  in  adaptive 
polarization  filter.  For  example,  it  is  easily  known  that  the 
estimation  precision  will  directly  limit  the  whole  filter 
performance  from  the  model  of  open  loop  cascading  system 
and  the  filter  performance  formula.  So  optimizing  the 
estimator  should  be  important  in  designing  the  filter. 

So  far,  adaptive  polarization  filters  are  mainly  used  to 
reject  jamming  with  high  polarization  degree.  Routine 
adaptive  polarization  filters  aren’t  satisfying  in  suppress 
jamming  partially  polarized.  This  is  a  desiderate  problem  in 
polarization  technology,  and  it  pose  new  demand  in 
evaluating  the  performance  of  adaptive  polarization  filters. 
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